Causality in Condensates: Grey Solitons as Remnants of BEC Formation. 
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Symmetry breaking during phase transitions can lead to the formation of topological defects 
(such as vortex lines in superfluids). However, the usually studied EEC's have the shape of a cigar, 
a geometry that impedes vortex formation, survival, and detection. I show that, in elongated traps, 
one can expect the formation of grey solitons (long-lived, non-topological "phase defects" ) as a result 
of the same mechanism. Their number will rise approximately in proportion to the transition rate. 
This steep rise is due to the increasing size of the region of the BEC cigar where the phase of the 
condensate wavefunction is chosen locally (rather than passed on from the already formed BEC). 
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Phase transitions are usually studied as equilibrium 
phenomena. However, as a consequence of the critical 
slowing down, second order phase transitions depart from 
equilibrium near the critical point, where the new bro- 
ken symmetry phase is chosen. Hence, that choice must 
be made locally, within regions that can dynamically 
"agree" on how to break symmetry. Cosmology offers 
a well-known example: As pointed out by Kibble [1 , rel- 
ativistic causality alone limits the size of domains over 
which symmetry breaking can be coordinated. As a con- 
sequence, topological defects such as monopoles, cosmic 
strings, and domain walls can form. 

In laboratory phase transitions relativistic causality 
does not provide useful estimates of the domain size with 
the approximately uniform new phase, and, hence, does 
not lead to predictions of defect density. One can, how- 
ever, estimate the domain size by appealing to univer- 
sality of second order phase transitions Symmetry 
breaking is coordinated by the dynamics of the order pa- 
rameter. In the vicinity of second order transitions criti- 
cal slowing down implies that the relaxation time (which 
determines the reflexes of the system) and the healing 
length (which sets the scale on which its order parameter 
"heals", i.e., returns to its equilibrium value) diverge as: 
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Above, To and depend on the microphysics, while the 
critical exponents v and z define the universality class of 
the transition, and e is the relative temperature 
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with Tc the critical temperature. 

Taking the ratio of ^ and r one obtains speed of sound: 
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This is the speed of perturbations of the order parameter. 
The resulting sonic horizon plays a key role. 

Divergence of the healing length was recently observed 
in measurements of phase coherence above the BEC crit- 
ical point [51. In effect, the experiment of Esslinger et al. 



demonstrated that the phase of the condensate wavefunc- 
tion is becoming coherent over distances that increase as 
the critical point is approached from above, as expected 
from Eq. ^ . If the critical region was traversed infinites- 
imally slowly, all of the newly created BEC would have 
a single coherent phase. However, when the transition 
is accomplished at a finite rate, critical slowing down, 
Eq. ([T|), intervenes: As its refiexes deteriorate, the phase 
of the order parameter cannot establish coherence over 
scales larger than the sonic horizon. 

In the usual discussions of topological defect formation 
[21 HI one first calculates the instant t at which the sys- 
tem ceases to follow the externally imposed variation of 
its parameters by comparing the timescale e/e at which 
relative temperature changes to the relaxation time: 



T{i) = e{i)/e{i) 
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To obtain t we need the dependence of e on t. We assume 
that it is linear, parametrized by quench time tq. 
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The system adjusts its state adiabatically as long as the 
imposed rate of change is slow compared to its reflexes 
given by the inverse of r, Eq. ([T]). The transition from 
the adiabatic to impulse behavior happens at t given by 



Eq. (5), i.e. tq\ 
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t . So, the order parameter 



"freezes" when the relaxation time and t coincide; 
t=(ror-)T^=f. 
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The order parameter will resume evolution only t after 
critical point is passed. The scale of the fluctuations (re- 
ported in Ref . [3] ) that seed structures (such as topologi- 
cal defects) in the broken symmetry BEC phase [2l[4l|5l|6] 
is thus established at t, i.e., at the relative temperature: 



The scale given by the corresponding healing length 
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determines the density of defects. The phase of the newly 
formed BEC wavefunction will be coherent on scales ~ ^. 
Therefore, one expects a defect fragment (e.g., one sec- 
tion of a vortex line) per ^-sized domain [T]. In a homo- 
geneous 3D quench this leads to vortex line density of 
~ [21 mil], in accord with most of the experimental 
evidence [B], including EEC's [HIS]- It is confirmed and 
refined by numerics jlOj . which also indicate that there 
is typically less than one defect fragment per ^-sized do- 
main: Rather, a defect fragment in /^-sized region, where 
rg-independent / is set by microphysics of the transition, 
is typical. The factor / can be greater than 1, and / ~ 10 
are common [10] . The density of defects created by phase 
transitions is the best known (but not the only) predic- 
tion of this "Kibble- Zurek mechanism" (or "KZM"). 

In the inhomogeneous case (e.g., effectively ID trap) 
situation is different: The gas density (and, hence, local 
critical temperature Tq) depends on location. Thus, even 
when T drops uniformly due to evaporative cooling, the 
gas will reach local critical temperature Tc (f) at different 
instants: e{fp,tp) — defines the front of the transition 
rpitp) as it spreads through the cigar. So the critical 
front will appear a,i tp that depends on location rp. 

Before the evaporative cooling the local density is [TT] : 



p{x) = poexp(-/3y(r)) 



(10) 



Above V{r) is (typically, harmonic) trap potential and 
(3 — l/ksT. Einstein's condition for BEC formation in- 
volves density and de Broglie wavelength, /3A^3(Tc) ~ 
2.61. In elongated traps one can in effect eliminate trans- 
verse dimensions [12j . This implies a local Tc {x) : 
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where m is the mass of bosons, while h and fcs are Planck 
and Boltzmann constants. In other words, when any- 
where in a large effectively ID harmonic trap tempera- 
ture falls below local Tc{x) in a region large compared 
to the healing length, the condensate will begin to form. 
We assume that cooling decreases T uniformly, so that; 



T{t) = Tc(0) 1 



t 
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everywhere in the trap. Therefore, front coordinates xp 
and tp are related by the equation e{xp,tp) = 0, or: 
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So the condensate can form first in a healing length size 
domain near a; = 0, where the potential is deepest. In an 
infinitesimally slow quench that initial seed would grow 
to occupy the whole trap. But this cannot happen when 
quench is so fast that regions far away from the center 



quickly attain temperatures far below the local Tc'{x): 
They will begin to form BEC independently, from local 
seeds, and with locally selected phases. 

The phase of the newly formed BEC wavefunction can 
be then either communicated along x, or selected at dif- 
ferent points of the trap independently (as would be the 
case in a homogeneous quench) . What actually happens 
is decided by causality, and depends on the comparison of 
the causal horizon defined by the relevant sound velocity, 
Eq. Q, and the velocity of the front: 



Vp 



dxf 



dtp 



Tc jO) 



dTcix) 



dxf 



(14) 



The speed of the front is infinite at the center of the 
trap where V{x) has its minimum, and drops with the 
inverse of the gradient of the critical temperature. The 
perturbations travel distance ~ ^ over time t. So, the 
relevant speed of sound corresponds to the freezeout e: 



(15) 



The role of v and its sonic horizon emerged in discussions 
of vortex formation in '^He superfluid. These experiments 
start with a cigar-shaped bubble heated above the crit- 
ical point [13] which quickly cools to the temperature 
of the surrounding '^He superfluid. One might have ex- 
pected that superfluid on the outside of the bubble will 
impose (uniform) phase of its wavefunction on the cool- 
ing "cigar". That this need not happen was noted in 
Ref. [14]: When the front T{x) = Tc{x) spreads faster 
than V, the phase of the newly formed condensate is cho- 
sen locally. Subsequent studies [TS] confirmed that when 
the front velocity vp exceeds u, symmetry breaking hap- 
pens as in a homogeneous transition, and defects appear 
with density inferred from ^. However, when v > vp, pre- 
existing condensate propagates its phase into the newly 
forming regions, and topological defects do not from. 

In the quasi- ID traps one does not expect to see vor- 
tices as, at formation, their ^-sized cores barely fit inside 
the cigar, so, even if they form, they can easily escape. 
Vortices do form in quasi-2D pancake traps [3 E] . But in 
the effectively ID geometry there is a stable defect related 
to phase nonuniformity - the grey soliton |16| . It corre- 
sponds to a solution of the Gross-Pitaevski equation, and 
describes a localized (healing length scale) nonuniformity 
of BEC phase, and a corresponding depletion of conden- 
sate density. The solitons are not topological: The phase 
change across the soliton can be arbitrary but, far (~ ^) 
away, it asymptotes to a constant value. Its change by 
TT yields a dark soliton, which causes complete depletion 
of the BEC density at its center. Dark solitons are sta- 
tionary, but grey solitons, with phase change less then tt, 
and a smaller depletion of central density (hence "grey" 
in their name), move along BEC cigar with velocities set 
by the local density depletion. When they arrive at the 
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point where BEC density is lower, they become (locally) 
"dark" , stop, and are reflected. Grey solitons were seen 
oscillating in this manner along BEC cigars [T7] . 

We can expect that non-uniformities of phase left by 
the BEC formation will give rise to grey solitons. Using 
KZM we can estimate density of phase jumps caused by 
the BEC formation. Thus, we can also estimate density 
of grey solitons in a homogeneous region, and their total 
number left in the trap by the phase transition into BEC. 

To this end, we compute local "freezeout" values of ^, 
f, and V using the local rate of change of e: 



de{x, t) . 
dt ' 
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This defines effective local quench time 

Tc{x) 
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which in turn yields local relative temperature: 

t — tp{x) 
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Now one can proceed as usual and compute local t: 

1 

=(ror5^(x))^ (19) 
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Note that this gives the time interval to the instant 
tp(x) at which the critical point is reached at the location 



X, and Eq. ( 16 1 is satisfied. This corresponds to the local 
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We have now all the ingredients to calculate the local 
frozen out healing length: 
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Local velocity at the freezeout is then a function of x: 
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These estimates are essentially the same as for the ho- 
mogeneous case: Key modification enters through the 
locally defined tq(x), Eq. (17 1. 

These predictions apply to the central part of the cigar 
where the critical front spreads faster than v - than the 
velocity of the perturbations of the order parameter. The 
region where the above quasi-homogeneous quench pre- 
dictions are accurate must therefore satisfy vp > v.^.. In 
view of our above discussion this leads to: 
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FIG. 1: Formation of grey solitons in a cigar-shaped Bose- 
Einstem Condensate, (a) Isodensity contour in the trapped 
gas. As evaporative cooling proceeds, critical temperature is 
first reached in the center of the trap. That is where the con- 
densate will form first. When cooling is sufliciently slow, this 
initial seed grows and imposes its selection of the condensate 
wavefunction phase on the whole cigar, and no grey solitons 
are created by the quench, (b) As is seen in the schematic 
color plot of the wavefunction phase in a cross section of a 
BEC cigar, the situation changes when BEC phase front - the 
location where the decreasing temperature is instantaneously 
equal to the local critical temperature (set by the local den- 
sity via Einsten's condition) - moves faster than the velocity 
i} with which perturbations of the emerging order parameter 
can spread. In this case regions of size ^, the relevant healing 
length, select phase of the BEC wavefunction independently. 
The front velocity vf is infinite at j: = 0, but falls rapidly with 
the distance x from the center. Condensate phase will be se- 
lected randomly by the symmetry breaking process in regions 
where v < vf- Such random phase distribution provides seeds 
for grey solitons. Phase front moves much less rapidly in the 
narrow direction of the cigar, so phases selected near the axis 
spread sideways, resulting in phase stripe pattern seen above 
in the schematic view of the BEC "cigar" . 



When V{x) = Tc{x), Eq. ([ll]), is a Gaussian: 

Tc(a;)=rc(0)e-^'/2^' (24) 
where A^^ = |/3rna;^, and we ignored variations per- 



penicular to the long axis. Inequality (23 1 leads to: 
Co tq' 



(25) 



This inequality determines size of the section [—X, X] of 
the cigar where vp > v, and the motion of the critical 
point is supersonic. There the quench is effectively ho- 
mogeneous, and defects (including solitons) will appear 
with separations given by the local ^ (see Fig. 1). 
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The equation for X is simple, but it is transcendental. 
We focus on the case where X < A. Then the exponent 



in Eq. (25 1 can be expanded, which leads to: 
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This estimate of X is valid for slow quenches, i.e. it 
breaks down when >i but holds when: 
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We assume that this is indeed the case. This focus on 
slow quenches is anyway prudent: Our discussion as- 
sumes that, outside of the freezeout interval, order pa- 
rameter is at or near the equilibrium set by the relative 
temperature e(t). Very rapid quenches could strain this 
assumption (although it is unlikely they could be imple- 
mented using evaporative cooling, as collisions that con- 
trol its rate also assure evolution of the order parameter) . 

Equation ( 27 1 yields simple scaling for the total num- 
ber of solitons. Note that above we have set = ^dB- 
de Broglie wavelength at the critical temperature [7 . We 
are now ready to estimate the total number of grey soli- 
tons. We obtain it by multiplying the size of the quasi- 
homogeneous quench region by the expected density of 
phase changes. This yields: 
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The surprise is that scaling of the number of solitons with 
the quench timescale tq is so steep. For example, for the 
plausible values — ^ and z — ^ we predict = |, 

while for mean field = 5 and z = 2 the exponent 
= 1. So the number of grey solitons is expected 
to be approximately proportional to the quench rate. 

Several aspects of the above prediction deserve com- 
ment. To begin, note that we have ignored all the aspects 
of the process that cannot be deduced from the universal- 
ity class. They will influence size of /. Here we include 
issues such as how dark a grey soliton must be to count 
as a soliton, and other matters relevant for experiments. 
For instance, it is known that solitons - while they are 
long-lived - do not live forever. Therefore, the number 
of solitons will depend on their survival rates. 

The calculation above also addresses the question of 
when the quench can produce a uniform BEG. This will 
happen when £^ ^ X, for quenches so slow that they 
produce N 1 solitons in a trap. There is also an 
opposite limit of very fast quenches. We shall not address 
it here as it is cumbersome (e.g., transcendental Eq. (25) 



cannot be approximated in a way that yields a simple 
result). Moreover, in order to reach it in experiments one 
would need to drop temperature very quickly throughout 
the trap, and far below Tc at the center of the trap. 



For such rapid quenches linear approximation e = t/rg 
is likely to break down, leading to further cumbersome 
but trivial complications. This last comment brings one 
more remark: In most experimental settings T(t) will fall 
below Tc at the center of the trap, but this may be above 
Tc{x) at some sufficiently large x. Our analysis applies 
as long as that x lies outside the central interval of 2X, 
or, more precisely, as long as the quench can be well 



approximated by linear relations (e.g., Eq. (12 1 inside it 



We discussed formation of grey solitons in elongated 
traps. The experiment aimed at detecting such non- 
topological remnants of a BEG phase transition should 
be easier than experiments [9] that study formation of 
vortex lines (which require more that a quasi- ID geome- 
try) . It should allow one to probe the connection between 
causality and symmetry breaking and test scalings pre- 
dicted by the Kibble-Zurek mechanism. 
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